Maximum gravitational recoil 
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Recent calculations of gravitational radiation recoil generated during black-hole binary mergers 
have reopened the possibility that a merged binary can be ejected even from the nucleus of a massive 
host galaxy. Here we report the first systematic study of gravitational recoil of equal-mass binaries 
with equal, but counter-aligned, spins parallel to the orbital plane. Such an orientation of the spins 
is expected to maximize the recoil. We find that recoil velocity (which is perpendicular to the orbital 
plane) varies sinusoidally with the angle that the initial spin directions make with the initial linear 
momenta of each hole and scales up to a maximum of ~ 4000 kms~^ for maximally-rotating holes. 
Our results show that the amplitude of the recoil velocity can depend sensitively on spin orientations 
of the black holes prior to merger. 

PACS numbers: 04.25.Dm, 04.25. Nx, 04.30.Db, 04.70. Bw 



Introduction: Generic black-hole-binary mergers will 
display a rich spectrum of gravitational effects in the 
last few orbits prior to the formation of the single rotat- 
ing remnant hole. These effects include spin and orbital 
plane precession, radiation of mass, linear and angular 
momentum, as well as spin-flips of the remnant horizon. 
Thanks to recent breakthroughs in the full non-linear nu- 
merical evolution of black-hole-binary spacetimes P, 0, 
3], it is now possible to accurately simulate the merger 
process and examine these effects in this highly non-linear 

regime 0, i S 0, i i E Ei E [11 Q M, El, [13, El. 

Black-hole binaries will radiate between 2% and 8% of 
their total mass and up to 40% of their angular mo- 
menta, depending on the magnitude and direction of 
the spin components, during the merger [i, [3, [i- In 
addition, the radiation of net linear momentum by a 
black-hole bi nary leads to the recoil of the final remnant 
hole [ii, [13, [m, m m, [H, [H, [H, [13, [13 . This phe- 
nomenon can lead to astrophysically important effects 

A non-spinning black-hole binary will emit net linear 
momentum parallel to its orbital plane if the individual 
holes have unequal masses. However, the maximum recoil 
in this case (which occurs when the mass ratio is q ^ 
0.36) is a relatively small '-^ 175 kms""'^ [23 |. 

The first generic simulation of black-hole binaries with 
unequal masses and spins was reported in [2J]. These 
black holes displayed spin precession and spin flips, and 
for the first time, recoil velocities over 400km s^^, mostly 
along the orbital angular momentum direction, ft was 
thus found that the unequal spin components to the re- 
coil velocity can be much larger than those due to unequal 
masses, and that comparable mass, maximally spinning 
holes with spins in the orbital plane and counter-aligned 



would lead to the maximum possible recoil. This max- 
imum recoil will be normal to the orbital plane. Brief 
studies of this configuration (with a/m between 0.5 and 



0.8) were performed in [2J, |26|. In this letter we report 



on the first systematic study of such configurations. Con- 
sistent and independent recoil velocity calculations have 
also been obtained for equal-mass binaries with spinning 
black holes that have spins aligned/counter-aligned with 
the orbital angular momentum ^23*, '25'| . Recoils from the 
merger of non-precessing black-hole binaries have been 
modeled in [3. 



In [2j| we introduced the following heuristic model for 
the gravitational recoil of a merging binary. 
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where ^ = 1.2 x 10^ kms-\ B = -0.93, H = (7.3 ± 
0.3) X 10"^ kms~^, di — Si/ra^, Si and mi are the spin 
and mass of hole i, q \s the mass ratio of the smaller 
to larger mass hole, the index _L and || refer to perpen- 
dicular and parallel to the orbital angular momentum 
respectively, ei, 62 are orthogonal unit vectors in the or- 
bital plane, and ^ measures the angle between the "un- 
equal mass" and "spin" contributions to the recoil veloc- 
ity in the orbital plane (see [2^ for a similar empirical 
formula). The angle 8 was defined as the angle between 
the in-plane component of A = m{S2/m2 — Si/mi) and 
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the infall direction at merger. We determine below that 
K = (6.0 ± 0.1) X 10^ kms"^ We note that the maxi- 
mum of the recoil velocity shifts toward equal-mass bina- 
ries when spin is present. For example, in the case where 
3.2 = —di = a the maximum recoil occurs for q — 1 both 
when = for acos(^) < 0.0 and when a" = for 
acos(e- Oo) > 0.07675. 

Current techniques are not accurate enough to mea- 
sure the spin directions of the individual holes at merger. 
Instead, we focus on the angle -d between the initial A 
(which, for our binaries, is parallel to the individual spins 
and to the orbital plane) and the initial linear (orbital) 
momenta of the holes. We test the dependence of the re- 
coil on 1? by varying the initial spin directions while keep- 
ing the initial puncture positions and momenta fixed. In 
addition, we choose configurations that suppress Vm and 
v± and maximize W|| . 



TABLE I: Initial data parameters. The punctures are located 
along the a;-axis at x/M — ±3.28413 with momenta P — 
±(0,0.13355,0), spins 5* — ±{Sx, Sy,0), and puncture mass 
parameters rrip/M = 0.430213. In all cases the specific-spin 
of the two holes is a/m — 0.5. 



Techniques: We use the puncture approach [3l| along 
with the TwoPuNCTURES dil thorn to compute ini- 
tial data. In all cases below, we evolve data containing 
only two punctures with equal puncture mass parame- 
ters, which we denote by rUp. We evolve these black-hole- 
binary data-sets using the LazEv [s^ implementation of 
the 'moving puncture approach' which was independently 
proposed in 0, Q . In our version of the moving puncture 
approach we replace the BSSN [s^, IHj Isll confor- 
mal exponent </>, which has logarithmic singularities at 
the punctures, with the initially field x = exp(— 4(/)). 
This new variable, along with the other BSSN variables, 
will remain finite provided that one uses a suitable choice 
for the gauge. An alternative approach uses standard fi- 
nite differencing of ^ We use the Carpet [sl] mesh 
refinement driver to provide a 'moving boxes' style mesh 
refinement. In this approach refined grids of fixed size 
are arranged about the coordinate centers of both holes. 
The Carpet code then moves these fine grids about the 
computational domain by following the trajectories of the 
two black holes. 

We obtain accurate, convergent waveforms and horizon 
parameters by evolving this system in conjunction with a 
modified 1+lqg lapse and a modified Gamma-driver shift 
condition 0, [s^, and an initial lapse a ~ 4'b'l (here 
ipBL = 1 + mp/(2ri) + TOp/(2r2), where is the coor- 
dinate distance to puncture i). The lapse and shift are 
evolved with (5* - /3*5,)a = -2aK, 5t/3" = B", and 
dtB°^ = 3/49tr" — r/i3°. These gauge conditions require 
careful treatment of Xj the inverse of the three-metric 
conformal factor, near the puncture in order for the sys - 
tem to remain stable 0, 0, [13] ■ As was shown in Ref . [s^ , 
this choice of gauge leads to a strongly hyperbolic evolu- 
tion system provided that the shift does not become too 
large. 

In Ref. [31 we presented convergence and consistency 
tests for our mesh-refinement code and showed that our 
code produces fourth-order accurate waveforms for spin- 
ning binaries with spin magnitudes equal to the spins 
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Sx 


Sy 


Madm/M 


SP2 








0.12871 


1.00001 


SPA 


-tt/4 


0.091013 


0.091013 


0.99980 


SPB 


Tv/2 


-0.12871 





0.99958 


SPC 


TV 





-0.12871 


1.00001 


SPD 


7r/2± 0.162 


-.12703 


-0.020760 


0.99959 


SPE 


-7r/2 


0.12871 





0.99958 



used in the present configurations and central resolutions 
of Af/32, M/40, M/52. 

Results: We evolved the configurations given in Ta- 
ble U with 9 levels of refinement and a finest resolution 
oi h = M/40. The outer boundaries were located at 
320M. We measure the gravitational recoil by analyzing 
the {£, m) modes of V'4 (^ < 4) as measured by observers 
at r 25M, 30M, 35M, 40M and extrapolating to in- 
finity. We take the error in our measured recoil to be 
the differences between a linear and quadratic extrapo- 
lation (in l/r) of these measurements. We removed the 
contribution of the initial (non-physical) radiation burst 
(which is typically 20km s""'^) from the computed re- 
coil. There are additional (small) errors due to our not 
including the initial non-zero recoil of the system as well 
as finite difference errors. Note that these configurations 
all have 7r-rotation symmetry, and consequently, if punc- 
ture 1 is located at {xp,yp, Zp) then puncture 2 will be 
located at (—Xp,—yp,Zp) (note the sign of the z coor- 
dinate). Thus these binaries do not undergo a typical 
orbital precession, rather the orbital plane itself moves 
up and down the z-axis. 

We obtained the momentum and puncture position ini- 
tial data parameters for the SP2 configuration using the 
3PN equations for a quasi-circular binary with period 
Muj ~ 0.0500 and the given spin. We determined the 
puncture mass parameter by requiring that the ADM 
mass be 1. We then rotated the spin, keeping its magni- 
tude constant, to obtain the parameters for the remaining 
configurations. 

Table |TT] and Fig. [1] give the recoil velocity for each 
configuration. A linear least-squares fit for all configura- 
tions yields = 1875.87 cos(i? - 0.183978). Note the 
very similar values for the radiated energy and angu- 
lar momenta. The fact that these values are identical 
to within 3% in the radiated energy and to within the 
errors in the calculation for the radiated angular mo- 
mentum, indicates that the binaries have very similar 
orbital dynamics. This expectation is further supported 
by the puncture trajectories in the xy-plane (see Figs. [21 
131 and|5]) which show essentially identical orbital trajec- 
tories and with an identical number of orbits prior to 
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TABLE II: The radiated energy and angular momentum and 
recoil velocity iiy for the configurations in Table |T] including 
predicted velocities based a least-squares fit of all configura- 
tions. 



Config nil (kms^^) v^^ (fit) B^ad/M 



SP2 
SPA 
SPB 
SPC 
SPD 
SPE 



1833 ± 30 
1093 ± 10 
352 ± 10 
-1834 ±30 
47± 10 
-351 ± 10 



1844 

1061 

343 

-1844 

41 

-343 



(3.63 ±0.01)% 
(3.53 ±0.01)% 
(3.57 ±0.01)% 
(3.63 ±0.01)% 
(3.55 ± 0.02)% 
(3.57 ± 0.02)% 



Jrad/M" 

0.248 ± 0.003 
0.244 ± 0.003 
0.246 ± 0.004 
0.249 ± 0.003 
0.245 ± 0.005 
0.246 ± 0.003 




FIG. 1: The recoil velocity versus angle between the initial 
individual momenta and spins and a least-squares fit. Note 
that the i? = ±7r are the same SPC configuration. 



merger. Thus we expect that a rotation of the initial 
spins will lead to an essentially identical merger but with 
the spins at merger rotated by that angle. Finally, we 
note that the PN spin-precession equations [iOl imply 
that the in-plane spin precession frequency is indepen- 
dent of spin orientation for our configurations to 1.5 PN 
order. Hence our fit of uy to cos(z9 — ^q) indicates that 
varies as cos(8 — Qq) as predicted by our empirical 
formula ([T]). 

Discussion: In an earlier paper 2J| we reported the 



first results from evolutions of a generic black-hole bi- 
nary, i.e. a binary containing unequal-mass (2:1) black 
holes with misaligned spins. These results suggested that 
the recoil velocities of rapidly-rotating black holes would 
be dominated by the contribution from the spins. While 
the configuration evaluated in that paper was not se- 
lected in order to maximize the recoil, the results were 
used to estimate the maximum recoil velocity due to 
spin, based on an empirical formula, Eq. ([T]). In this 
Letter, we have confirmed our previous estimates with a 
set of new numerical simulations of binaries having spins 
of equal magnitude but counter-aligned, and parallel to 
the orbital plane. We found that these configurations 
maximize the W|| term in Eq. ([T]) while setting the re- 
maining terms to zero. We confirmed that U|| varies as 



-1.5 - 



-3.5 




x/M 



FIG. 2: The projection of the puncture trajectories (only 1 
shown per configuration) for the 6 configurations. The orbital 
dynamics of the binaries are not significantly affected by the 
change in spin directions. 




FIG. 3: The z-component of the punctures trajectories (only 1 
shown per configuration) versus time for the 6 configurations 
showing the dependence of the orbital plane 'precession' and 
remnant recoil on the angle of rotation. 



K cos{{} — i?o), where d measures the angle between the 
initial spin and linear momentum vectors. Based on the 
fit U|| = (1876 ± 30) cos(z?- 0.183978) we determined that 
K = (6.0 ±0.1) X 10*kms~^. Since the magnitude of 
the recoil predicted by Eq. ([1]) is proportional to the di- 
mensionless spins di , our results predict maximum recoil 
velocities of ^ 4000 kms~^ in the case of maximally- 
spinning holes with counter-aligned spins. 

A post-merger recoil velocity of ~ 4000 kms"-'^ is large 
enough to eject a black hole from the center of even the 
most massive elliptical galaxies [s^l- Hence, our results 
strengthen the conclusion, already reached in several re- 
cent papers 24, 2^, 28 1 that radiation recoil is capable of 
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FIG. 4: The three-dimensional trajectories of the punctures 
showing the orbital precession and the final recoil for the SP2 
configuration. Note that the scale of the z-axis is 1/10 that 
of the X and y axes. 



completely removing supermassive black holes (SMBHs) 
from their host galaxies. Computing the probability of 
such an extraordinary event will require a more extensive 
set of numerical simulations that characterize the depen- 
dence of V^ocoii on spin direction for generic binaries, with 
arbitrary spin orientations and mass ratios. Here, we 
note the strong predicted dependence of Kocoii ^ 9^ on 
mass ratio which implies a "maximum" recoil velocity of 
~ 0.3^ X 4000 kms~^ « 400 kms~^ even for a "major 
merger" with m2/mi « 1/3. In addition, the root-mean- 
square recoil velocity for randomly oriented spins in the 
plane is reduced by an additional factor of ^/2. Our re- 
sults are therefore not inconsistent with the observed fact 
that SMBHs are apparently ubiquitous components of lu- 
minous galaxies. A galaxy in which the SMBH has been 
permanently removed will appear similar to a more or- 
dinary galaxy except for a larger stellar core, or "mass 
deficit," due to heat input from the ejected hole. In fact, 
observed mass deficits are sometimes ^ a few times larger 
than predicted on the basis of binary SMBH models 4l| . 
The detection (or non-detection) of a SMBH in such a 
galaxy via stellar absorption line features will be difficult 
however due to its low central surface brightness. 
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